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ABSTRACT 

Further  analytical  investigations  are  made  into 
the  damping  capability  and  determination  of  natural 
frequencies  of  laminated  beams,  consisting  of  elastic- 
viscoelastic-elastic  layers,  as  a  means  for  reducing 
the  vibratory  energy  transmitted  through  machine  founda¬ 
tion  supports  in  naval  vessels . 

An  exact  analytical  solution  is  obtained  for  de¬ 
termining  the  natural  frequencies  of  simply-supported 
sandwich  beams  having  no  rivets  at  the  ends0  Three 
possible  modes  of  vibration  are  shown  to  exist,,  The 
case  of  the  simply-supported  sandwich  beam  having  rivets 
at  each  end  is  considered  and  the  equations  reduced  to 
the  solution  of  12  x  12  determinant  for  calculation  on 
a  digital  computer. 

An  approximate  method  is  suggested  for  determining 
the  natural  frequencies  of  sandwich  beams  having  any 
end  conditions.  The  procedure  is  simple  to  use  and  is 
exact  for  simply-supported  beams. 

A  simple  but  approximate  expression  is  also  de¬ 
veloped  for  determining  the  composite  loss  factors  of 
sandwich  beams.  The  procedure  yields  good  engineering 
results. 
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Nomenclatures 

5-j  =  mid«plane  extension 
R  =  R-j  (Hig) 

3  =  material  loss  factor 
G-jb 

Ri =  agq- 

b  =  width  of  beam 

G  =  G^(Hie)  -  complex  shear  modulus 
G-j  =  real  part  of  complex  modulus 
2H-j  =  thickness  of  elastic  layer 
2H2  =  thickness  of  viscoelastic  layer 
K  =  EA 

E  -  elastic  modulus 
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Natural  Frequencies  and  Damping  Capability  of  Laminated  Beams 
Io  Synopsis 

loO  Introduction 

1.1  Overall  Program;  As  part  of  the  overall  Navy  Program  to 
reduce  noise  emanating  from  vessels,  a  "Structural  Damping 
Program5  is  in  progress  at  the  U,  S„  Navy  MEL,  The  aim 

of  this  program  is  to  investigate  methods  for  attenuating 
vibratory  energy  in  the  structure  between  the  machines  and 
the  hull. 

1.2  Specific  Program;  The  UoSoNoMoEoL0  has  taken  an  approach 
to  accomplish  the  above  goal  through  the  use  of  laminated 
material  as  structural  members,,  This  material  is  made  of 
steel  and  viscoelastic  layers,  so  that  structural  rigidity 
is  maintained  concurrent  with  damping  (or  dissipation  of 
vibratory  energy )a 

103  Previous  Results-  Tests  have  shown  that  reasonable  amounts 
of  damping  are  possible  by  using  laminated  material „  In 
order  to  better  understand,  use  and  design  with  the  material 
arialytica1  investigations  have  been  in  progress,,  The  re* 
suits  of  References  1,2,3  and  4  indicate  that  mathematical 
expressions  are  available  for  determining  the  composite 
loss-factor  of  a  laminated  beam  as  a  function  of  the  natural 
frequency,.  These  results  also  show  that  the  composite  loss- 
factor  versus  frequency  curve  is  independent  of  the  boundary 
conditions  and  solely  dependent  upon  the  physical  and  geo¬ 
metric  properties  of  the  cross-section  of  the  beam,  The 
natural  frequencies  themselves  are  dependent  upon  the  boun¬ 
dary  condi tionSo  In  reference  5  the  free  and  forced  vibra¬ 
tions  of  sandwich  beams  were  investigated,.  The  eigenvalues 
no»  s: »  s2j  t;.  and  t^  were  obtained  as  a  function  of  h;  and 
plotted.  The  mathematical  expression  for  y,  <}>»  M  and 
V  were  explicitely  written  in  terms  of  the  eigenvalues„  A 
generalized  plot  is  given  of  the  composite  loss  factor  versus 
hi  and  ($  from  which  the  composite  loss-factor  n  can  be 
readily  found.  Reference  6  is  a  short  survey  of  the  theory 
of  viscoelasticity  especially  as  it  pertains  to  this  pro¬ 
gram., 

2o0  Scope  of  this  Study;  Details  of  the  free  vibrations  of 

sandwich  beams  have  been  investigated.  The  "exact9  natural 
frequencies  have  been  obtained  for  a  simply-supported  sand¬ 
wich  beam.  In  view  of  the  complexity  of  the  "exact"  solution, 
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approximate  solutions  have  been  evolved  which  can  be  used 
to  calculate  the  natural  frequencies  of  sandwhich  beams 
having  any  end  conditions,,  These  approximate  procedures 
allow  one  to  design  beams  for  optimum  damping  and/or 
measure  the  properties  of  the  viscoelastic  material „ 

3.0  Results 

3.1  Natural  Frequencies  -  Exact:  The  two  cases  considered 
are  the  simply-supported  sandwich  beam  with  and  without 
rivets  at  each  end*  The  unriveted  case  yielded  the  analy¬ 
tical  result  that  a  natural  frequency  exists  when 

sin  no"2  R,,'2  L  =  0 


or 


n0'2  R'2  L  =  mr  . 

This  result  is  to  be  expected  and  is  the  result  for  a 
sinusoidal  deflection  y  =  y0  sin  „  An  unexpected  re¬ 
sult  is  that  natural  frequencies  exist  when 

2tJ  L  =  rnt 


and 


2t?  Rf  L  =  imr  . 

This  result  appears  to  be  new  and  indicates  that  natural 
frequencies  can  occur  in  simply-supported  sandwich  beams 
for  other  than  sinusoidal  mode  shapes .  The  t{  vs  hi  curves 
show  that  only  the  lower  values  of  n  yield  natural  fre¬ 
quencies  and  a  cut-off  value  of  n  exists  above  which  none 
of  this  type  of  mode  of  vibration  occurs.  The  t^  vs  h2 
curve  displays  a  similar  characteristic;  i.e.  only  the 
first  few  values  of  n  yield  natural  frequencies  and  a  cut¬ 
off  value  of  n  exists  above  which  no  natural  frequencies 
exist. 

The  case  having  riveted  ends  did  not  yield  to  a  closed- 
form  solution  and  required  the  calculation  of  a  12  by  12 
determinant.  The  necessary  elements  of  the  determinant 
are  presented  in  a  form  ready  for  use  in  a  digital  com¬ 
puter. 
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3,2  Approximate  Relations  for  detaining  the  Natural  Frequencies 
of  Sandwich  Beams , 

(a)  The  natural  frequency  is  written  in  the  form 
a 

w2  =  *■— —  [B^  +  Bo  +  S2a] 

*  pL2 

in  which 


a  is  a  constant  dependent  upon  end-conditions 
{see  page  35), 

p  is  the  mass/unit  length 

lL 

B-j  is  the  stiffness  of  the  t  elastic  layer., 

Ki  =  EAj 

Ai  is  the  cross-sectional  area  of  the  elastic  layer,, 

6  is  the  distance  between  the  neutral  axes  of 
layer  1  and  3, 

a  is  shear  parameter  given  by  equation  JEK'tf) 

For  the  case  in  which  e  <<1,  it  is  found  that 


2H.2H 


a  -  [ap(- 


2  c 

-)  £■+  2]  1 


3o3  A  simplified  method  for  determining  the  shear  moduli  and  loss- 
factor  of  a  viscoelastic  for  B.«l  is  proposed  in  section  III 
in  which  the  composite  loss  factor  n  and  natural  frequency 
are  measured  and  Gj  and  6  calculated.  This  procedure  may  be 
used  as  a  procurement  test-method  for  the  viscoelastic  material. 


3,4  Optimum  damping  for  a  sandwich  beam, 

(a)  Design  -  Given  E,  H2;  y»  «i  and  the  end  conditions,  one 


can  determine  the  optimum  length  of  the  beam  and 
elastic  layer  thickness  to  yield  maximum  damping 
for  a  given  viscoelastic  material* 

(b)  Design  of  the  shear  properties  of  the  viscoelastic 
may  be  performed  so  that  optimum  damping  may  be  G] 
obtained  over  a  wide  frequency  range*  By  letting  — 
be  constant  for  the  maximum  value  of  n»  maximum 
damping  is  attainable  over  the  frequency  range* 

Thus,  if  the  real  part  of  the  shear  modulus  is  speci¬ 
fically  designed  into  a  material  to  increase  wlt  then 
an  optimum  material  can  be  obtained*  In  all  cases  the 
tangent  loss  factor*  should  be  as  large  as  possible* 

3*5  An  approximate  equation  has  been  found  for  the  composite 
loss  factor  of  a  sandwich  beams  which  is  within  10%  of 
the  correct  value*  The  relation  is 


n  *  6  6  6  — '  n*™™”  —  r\  "  O 

( (2+*475  h:D)^+{26)  3 

3*6  It  is  further  shown  in  section  III,  that  the  natural  fre¬ 
quency  of  a  sandwich  beam  increases  with  damping*  This  is 
in  agreement  with  a  result  derived  in  reference  3* 

4*0  Continued  Investigations 

It  is  recommended  that  the  following  be  considered  for 
further  investigations  of  sandwich  beams* 

4*1  Computer  solutions  be  performed  to  find  the  natural  fre¬ 
quencies  for  the  following  end-condition 

a*  Simply-supported  riveted  ends 
b*  Free -free 
c*  Cantilever 
d*  Fixed-fixed 
e*  Fixed-pinned 

Correlation  with  the  approximate  methods  should  be  investigated* 

4*2  Consider  the  mode  shapes  and  meaning  of  the  natural  frequencies 
associated  with 


tjL  -  nu  and  t2L  =  mir* 

Consider  correlation  with  tests  if  necessary 


Of>  ^  oo 


rr**? rr-~- 


lie  Determination  of  the  natural  frequencies  of  sandwich  beams  -  exact 
solution* 

1*0  As  discussed  in  reference  5,  the  determination  of  the  natural 
frequencies  of  three  layer  laminated  beams  may  be  found  by 
satisfying  six  boundary  conditions,  i*e*,  the  usual  four  of 
deflection,,  and/or  slope,  and/or  moment  and/or  shear;  plus 
an  extensional  effect  at  each  end,  i*e*  zero  stress,  and/or 
zero  deflection  and/or  zero  shear* 

The  equations  for  extensional  effect 

5  *  0 

4"  =  0 

C*  ’  *  0 

and  for  'lateral*  effects 
y  =  0 
<p  -  0 

M  *  0 
V  -  0 

are  given  in  reference  5,  as  equations  II-1C,  1 1 -2,  1 1-3, 

11*5,  II  7,  II-8*  and  1 1  *10.,  Equations  for  M  and  V  are  re* 
written  in  the  appends  as  errata  to  1 1-8  and  11  =  10*  The 
general  case  of  different  elastic  layers  can  be  solved  on 
the  basis  of  'individual  cases*  In  general,  the  solution 
requires  the  evaluation  of  a  12  by  12  determinant  and  the 
results  cannot  be  eas*iy  generalized*  The  special  but  im¬ 
portant  results  of  the  sandwich  beam  can  be  generalized  and 
this  case  has  been  solved  herein  for  the  simply-supported 
beam  with  and  without  rivets  at  each  end*  The  simply-supported- 
no-rivets  case  requires  that 

M  ~  y  =  4  -  0, 

since  the  moment  defection  and  axial  stress  is  zero  at  each 
end*  These  six  boundary  conditions  yield  12  linear  homogeneous 
equations  By  selecting,  the  coordinate  system  at  one  end,  the 
twelve  equations  become  two  sets  of  six  equations  each  having 
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its  own  set  of  six  undetermined  coefficients.,  Since  one  set 
is  independent  of  the  length  of  the  beam  (the  only  remaining 
variable  in  the  determinant  of  coefficients)  then  this  set 
can  only  be  satisfied*  in  general*  if  the  six  undetermined 
coefficients  in  that  determinantal  set  are  all  zero*  The 
determinant  of  coefficients  of  the  second  set  of  six  equations 
can  be  made  equal  to  zero  for  three  cases 

sin  a^L  *  0 

sin  2t^L  «  0 


and 


sin  2t2L  '  0  . 

The  first  case  's  the  expected  one  in  which  the  mode  shape  is 
sinusoidal  (y  =  yQsin  tip.)  and  the  natural  frequencies  are 
obtained  when  L 


or 


„  - 

nr>  ~  VT”"' 


1 

fT 


where  n  is  an  integer  1,  2,  3*  etc,  For  a  given  Lr  R-|  and  n 
the  value  of  n0  may  be  computed.  With  this  value  of  n0 
and  a  value  of  8*  one  may  find  h]  from  the  curve  on  page  11-34 
of  reference  5,  These  natural  frequencies  correspond  to  the 
usual  relation 


n2  2 

n  7r 


El 


o 


The  second  and  third  cases  were  unexpected  results.  These 
imply  that  natural  frequencies  also  occur  for  a  simply-supported 
no-rivets  case  at  other  frequencies  and  mode-shapes  than  sinu¬ 
soidal,  The  mode  shapes  were  not  investigated  but  the  frequencies 
were.  For  a  natural  frequency  to  exist 


2t-j L  -  nn 
or 

2t2L  =  mu 


—  6  — 


where  m  and  n  are  integers,.  A  given  value  of  n  or  m  arid 
L  yields  a  value  for  t]  and  t2  which  can  be  used  in  curves 
11-37  and  11-38  of  reference  5.  Recalling  that 


and 


the  values  of  tj  and  t£  may  be  calculated  and  used  in  the 
curves  to  obtain  h-j  from  which  wj  may  be  obtained0 

It  is  seen  that  the  tf  type  of  vibration  could  yield  an  h-j 
for  a  given  L  and  n  =  1  and  for  n  =  2  it  is  possible  for 
lower  value  of  h-j  to  be  obtained.  At  some  value  of  n  and  8 
the  curves  indicate  a  cut-off  of  any  of  these  type  of  vibra¬ 
tions.  Thus  the  lower  mode  is  associated  with  a  high  natural 
frequency  of  vibration  in  this  kind  of  vibration.  This  is 
contrary  to  the  usual  vibration  phenomenon  and  merits  further 
academic  investigation.  The  tj>  type  of  vibration  can  yield 
a  value  of  h]  for  m  *  1  and  a  given  L.  For  higher  values  of 
m,  larger  values  of  h-j  may  be  found  until  a  cut-off  frequency 
is  reached.  This  too  merits  further  academic  study.  The  simply- 
supported  beam  with  riveted  ends  requires  that 

M  =  0 
y  =  0 


and 


Ct0  =  0 

at  each  end. 

Since  each  aquation  is  complex,  the  six  boundary  conditions 
yield  twelve  linear  homogeneous  algebraic  equations.  Unlike 
the  unriveted  case,  the  resulting  12  x  12  determinant  does 
not  simplify  by  factoring.  In  this  case,  the  origin  of  the 
coordinate  system  was  placed  at  the  center  of  the  beam  and 
the  equations  written  accordingly.  In  order  to  retain  some 
generalization  to  the  results,  a  factor  of  R-j  was  eliminated 
from  each  of  the  equations  by  using 
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1, 

and  using  F  =  RfL  as  a  parameter,  The  specific  equations  and  the 
144  elements  of  the  resulting  determinant  have  been  written  and  are 
to  be  used  to  find  the  natural  frequencies  of  this  type  of  beam. 
Computations  are  in  progress  on  the  IBM  360  computer  for  finding 
the  first  five  natural  frequencies,,  Using  the  results  of  reference  5 
for  n0»  s{,  S£»  t{  and  t^  and  values  of  equal  to  6,  20,  50,  150, 

500,  2000  and  5000  for  0  equal  to  .1  and  1,  values  of  F  will  be  found 
which  satisfy  the  12  x  12  determinant.  These  results  will  be  reported 
in  the  future  as  a  plot  F  vs  hi  for  the  first  five  modes  of  vibration. 


oU  Simply-supported  ■  3  layer  beam  »  ends  unconstrained 

The  boundary  conditions  involve  the  Moment,  the  deflection  y  and  the 
axial  strain,, 

Now 

H  =  K,$c}  +  (B,  +  B3)  f  -1 

and 


=  f  ‘ 


1 

W 


iM 


-2 


so  that 

M  ■  tv  +  <W  N  -  Air1)  «i" 

also 

=  0  @  x  =  0  &  x  =  L 

and  in  addition 


In  view  of  the  fact  that  must  be  zero  at  each  end,  then  this 
term. in  M  and  y  may  be  dropped  so  that  for  these  boundary  conditions, 
there  remains 


/? dx  =  0 

C1  =  0 


@  x  =  0 
and 

x  =  L. 


In  order  to  evaluate  the  above  expressions  the- values  of  the  coefficients 
S  •  • ,  T. U. .  and  /R .  -  dx  are  evaluated  at  x'=*0s  They  are; 

*  J  *  v  ♦  J  *  •  J 
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R31 

=  1 

R32 

=  0 

R41 

=  0 

R42 

=  0 

R51 

=  1 

R52 

=  0 

R61 

=  0 

R62 

=  0 

S31 

=  0 

S51 

=  0 

S32 

=  0 

S52 

=  0 

S41 

=  S1 

S61 

=  S2 

S42 

=  tl 

S62 

=  t2 

T31 

=  (s 

2  t  2, 

1  _tl  l 

T51 

*  (s22-t22) 

T32 

=  2s 

ih 

T52 

^S2^2 

T41 

=  0 

T61 

=  0 

T42 

=  0 

T62 

=  0 

U31 

=  0 

U51 

=  0 

U32 

=  0 

U52 

=  0 

U41 

=  sl‘ 

:s12-3t12] 

U61  ; 

=  s2[s22-3t22] 

U42 

■  tjOs/.t,2] 

U62  : 

t2[3s2  -t2  ] 

/R31 

dx  = 

0 

^R51 

dx  =  0 

/r32 

dx  = 

0 

J"R52 

dx  =  0 

fR„, 

dx  = 

S1 

fo 

-  c. 

*  41 

77777 

si  n 

JR61 

dx  2  2 

s2  z2 

[RAn 

dx  = 

ti 

fD 

to 

<4  w  —  ^ 

•'42 

£  2 
si  +  *i 

JR62 

s2  +  t 
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* “u1  nn  ijp^""p'  t'.vir 


(D  at  x  =  0 


cj  =  0 

so  that 


(1-1)  E^i a-j  +  E41s1  -  E42t1  +  Egls2  -  Eg2t2  =  0 


and 


(1-2)  E22a^  +  E4i t-j  +  E42s-j  +  Eg-jtg  +  Eg2s2  -  0 

(7)  at  x  =  0 

c'"  =o 

so  that 

(2-1)  -E2^a-j  +  E 4 *j s i ( s ^  “3t-|  )  -  E42 1 1 ( 3s -j  -t-j  )  +  ^g]S2[s2  "*3^2  ^ 


-  E62t2(3s2  -t2  )  -  0 


and 


(2-2)  ”^22^1  ^41^*1  ( 3s i  -t-j  )  +  E42s-j(s.j  -3t-|  )  +  ^62s2^2  ^ 


2  04.  2\  .  r  _  /_  2_  2| 


+  E61t2^3s22’t22^  ~  0 


(7)  at  x  =  0 
Jc dx  =  0 
so  that 


(3_1)  'E21  af  +  e41  7T777  +  e42  7T777  +  e61  7  '2'  7 T  ? 

1  s,  +  t,  s,  +  t,  s2  +  t2 


+  E  2 

o2  - 7  =  0 

s  C  +  t  “ 
s2  z2 
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X 


(3-2)  -E„  i-  -  E 


1 


r  S 1  r-  ^2 

22  a7  "  l41  7  2  77?  +  E42  c  2  ,  .  £  ’  E61  TTTTT 
1  s1  +  t,  s,  +  t]  s,  +  t. 


’2  2 


+  £  ^ 

+  o2  “TTTT 


=  o 


S2  +  l2 


These  six  equations  are  all  independent  of  the  length  L,  indicating 
that  the  only  non-trivia!  solution  is  for: 


E 

E 

E 

E 

E 

E 


21 

22 

41 

42 
61 
62 


0 

0 

0 

0 

0 

0 


@  x  =  L 


R 

R 

R 

R 

R 

R 

R 

R 


31 

32 

41 

42 

51 

52 
61 
62 


=  cos  t-jL  cosh  s^L 
-  sin  t-jL  sinh  s-jL 
*  sinh  s -j L  cos  t^L 
=  sin  t-jL  cosh  S-jL 
«  cos  t2L  cosh  SgL 
=  sin  t2L  sinh  s2L 
=  sinh  s2L  cos  t2L 
=  sin  t2L  cosh  s2L 


531  =  s1R41  "  tlR42 

532  =  S1R42  +  tlR41 

541  =  S1R31  '  tlR32 

542  *  S1R32  +  tlR31 

551  =  S2R61  "  t2R62 

552  =  s2R62  +  V'ei 

561  =  S2R51  “  t2R52 

562  =  S2R52  +  ^51 


12 


=  s1  [s J2-  St^]  R41  -  t1  [3S-J2  - 
=  s i  [ s i  -  3t i  }  +  l*i  [ 3s -j 

=  s -j  [ s -j  -  3t-j  ]  R^ -j  ~  ti  [  3s -j 

*  s1  [S^-  3t12]  R32  +  t1  [3s-,2  - 

s?[s22-  3t2  ]  Rgi  -  t2  [3s2  - 

s2[s2  -  3t2  ]  Rg2  +  t2  [3s2  - 

=  s2[s2  -  3t2  ]  Rgi  -  t2  [3s2 

=  S2[S2  -  3t2  ]  Rg2  +  t2  [3s2 


t]2]R42 

tl2]R41 

tl2]R32 

tl2% 

t22]R62 

^^61 

^^52 


Letting 


We  have 


pij  =  /Rij  dx 


si  ‘i 

P-al  =  5  7  R/il  +  5  T  »  -p 


31 ' 7. ^7  52  -  *62  -  *61 


7T7TT  R31  +  "  2, .  2  R32  ;  P62 

5 1  '  t  j  5 i  •  u*j 

51  n  n 

T~±~~£  R32  "  "  5  ~  T  R31 
Si +t,  iL  s, +y 

52  ^2 

— z — TT-Rr-.  +  — * — R«-o 


_  fc.ifc  W  I  .  t,.  t 

S2  n  S2  +t2 


hW Rsi  +  ^7 R'2 

s2  ^2  n 
2  2  R52  5  ?  R51 

S2  +t2  s2  +t2 
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!:==^ssse**>^™**i'u»' 


Considering  the  boundary  conditions  at  the  other  end, 

(D  at  x  =  L 
0 

(4-1)  -E„  a,  sin  a,L  ♦  E3,S3,  -E^  +  E5,S5,  -E^  -  0 
and 

(4-2)  -E12  a,  sin  a,L  +  E3,S32  +  E32S31  +  E51S52  ♦  E^,  -  0 

CD  at  x  =  L 

c'"  =0 
so  that 

(5-1)  E„  a,  3  sin  a,L  +  E3]U3,  -E^  +  E5)U5,  -E^  =  0 
and 

(5-2)  E12  a,  3  sin  a,L  +  E3)U32  +  ♦  E^U^  +  E52U51  =  0 

(S  at  x  *  L 

/  t  dx  *  0 
so  that 


sin  a,L  s,  t, 

(6-l)  En  {  r  )  +  E31[-  -2  +  2  R41  +  §  '  .  2  R42] 

i  s1  +  z}  s1  +  t] 


"e32cs  l  \  t  2  R42  "  7T77T  R41]  +  e51[7T777  R61  + 
S1  S1  z]  s2  +  t2 


"£52[  T77T  R62  "  2  7  7  2  R613 

s2  +  t2  s2  +  *2 
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sin  a-, L  s,  t, 

(6.2)  e12(  — )  +  e31[  2  2  R42  -  ~T~7  k41j 

*1  +  *1  S1  +  *i 


R,  1 


■E32[7T77T  R41  +  7T777  R42] 
S1  'l  S1  +  tl 

r62  ■  7/77? Rfil1 


"E52[  i  \  2  R61  +  77 1"  2  R62!  *  0 

s2  +  *2  s2  +  E2 


form 

(6-3) 


It  is  seen  that  when  the  six  above  equations  are  written  in  matrix 


a..  E.. 

ij  ik 


=  0 


that  the  first  two  columnsare  made  up  of  a  the  common  factor  sin  a-jL .  so 
that  the  matrix  equation  may  be  written  as 


(6-4) 


sin  aL  Ib.^lE.J  =  0 


in  which 
(6-5) 


The  matrix  ( b^ ^ J  is  the  same  as  |a^|  except  for  the  first  two  columns 


an =  bn sin  aiL 


ij 


and 


(6-6) 


ai2  =  bi2  sin  alL* 


Equation  ( 
(6-7) 

or 

(6-8) 


or 


(6-9) 


)  shows  that  a  solution  is  obtained  when 
sin  a^L  =  0 


a-jL  =  mr 


,  _  nTr 

li  -  r 
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Vi'WwjSiLIJUI.. 


“r*X£2! XT***  TV' 


where 


.  determinant  has 

•J 

the  following  elements 

’n  =  _ai 

b21  3 

0 

o 

II 

04 

b22  3 

-al 

'13  =  slR4l“tlR42 

b23  3 

slR42+tlR41 

'14  =  ”slR42”tlR41 

b24  3 

SlR41-tlR42 

'15  =  s2R61't2R62 

b25  3 

s2R62+t2R61 

'16  =  *s2R62*t2R61 

b26  = 

S2R61“t2R62 

'31  =  213 

b41  3 

0 

'32  3  0 

b42 

ai3 

33  3  ClR4l”C2R42 

b43 

C1R42+C2R41 

'34  =  “C1R42“C2R41 

b44 

C1R4TC2R42 

35  3  C3R61“C4R62 

b45 

C3R6l”C4R62 

36  3  “C3R62+C4R6! 

b46  3 

C3R62+C4R61 

^1  ~  ^1^1  *1  ^ 

and  C2  -  t-,  (3s-j2-t-j2) 

S2  ( )  and 


^4  -  ^2^S2  ”^2  ^ 


> 


J51 


i 

I- 


b52  =  0 


b53  =  dlR41+d2R42 
b54  "~dlR42+d2R41 
hb  =  d3R61+d4R62 
b56  =  “d3R62+d4R61 


where 


dl  = 


£  ? 
s,  «,2 


and 


b61 

bB?  *  r- 

b63  =  dlR42"d2R41 

b64  =  ”dlR4l”d2R42 

b65  =  "d3R62+d4R61 

b66  =  ”d3R6l”d4R62 
ti 

d2  =  2  2 
si  s 
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^P*P*J*^  V 


-  4  , 


•> 


on-  (6-4)  may 

be  written 

as  the 

sum  of  two 

matrices 

or 

b! .  E,  .  + 

ij  kj 

b!i  E, 

1 3  kj 

=  0, 

al 

0 

s1R41 

"s1R42 

s2R61 

“s2R62 

E11 

0 

-al 

tlR41 

-tlR42 

*2^1 

“t2R62 

E12 

al3 

0 

C1R41 

-c1R42 

C3R61 

"C3R62 

E13 

0 

al3 

C2R41 

“C2R42 

C3R61 

C3R62 

E14 

0 

0 

dlR41 

-dlR42 

d3R61 

“d3R62 

E15 

0 

0 

“d2R41 

~d2R42 

d4R61 

"d4R62 

E16 

0  ' 

0 

-tlR42 

“V41 

”t2R62 

_t2R61 

E11 

0 

0 

s1R42 

s1R41 

S2R62 

"S2R61 

E12 

0 

0 

wC2R42 

"C2R41 

“C4R62 

C4R61 

E13 

0 

0 

C1R42 

C1R41 

”C4R62 

C4R61 

E14 

1 

al 

0 

d2R42 

d2R41 

d4R62 

d4R61 

E15 

0 

1 

al 

dlR42 

~dl R41 

“d3R62 

”d3R61 

E16 

=  0  (6*10) 


-  17  - 


It  is  seen  that  the  above  equation  may  be  written  in  the  form 


Tlvis  may  be  interpreted  in  the  fpl lowing  manner.  Given  a  beam  in  which  R-,  and  L 

G,n 

are  known  then  assume  n  =  l*a  value  for  tj  is  calculated,  (recall  R-|  =  ^|-y-) . 
This  value  of  t^j  is  used  to  find  h-j  for  a  given  6  using  the  chart  on  pg  II  -  37 
ref  (  5).  Knowing  h-j,  the  corresponding  composite  loss-factor  n  may  be  found 
from  pg  II  -  33  of  the  same  reference  and  the  natural  frequency  w-j  may  be  cal¬ 
culated  using 

hl  ;  48  h22  “I2  • 

G1 

Integer  values  of  n  (2,  3,  4  etc)  may  be  substituted  to  obtain  larger  values 
of  t|  and  h  /  use  of. the  charts  obtain  the  corresponding  natural  frequency 
and- composite  loss-factor.  For  a  given  s  the  t-j  vs  h-j  curves  indicate  lower 
natural  frequencies  for  increasing  n  with  a  cut-off  frequency  for  which  no 
natural  frequencies  occur  above  a  given  value  of  n. 

A  similar  procedure  may  be  followed  to  obtain  the  natural  frequencies  and 
associated  .composite  loss-factor  for  t£  using  the  curves  of  t£  vs  h^  on  pg  II  - 
38.  In  this  case  the  natural  frequencies  increase  with  m  and  reach  a  cut-off 
frequency  for  a  given  value  of  3. 
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3.0  Simply-supported-3  layer  beam-ends  constrained  (riveted) 

The  boundary  conditions  at  both  ends  are 

M  =  0 
y  =  0 

and  5“  =0,  (ends  riveted). 

The  moment  and  deflection  equations  do  not  simplify  as  for  the 
previous  case  with  no  rivets.  The  moment  equation  of  reference  & 
is  corrected  and  shown  as  errata  at  the  end  of  this  report. 

For  ease  of  writing^the  three  equations  are  written  as 

M  ■  +  “ 0 

y  *  ZZtdjjXj  +  l9jixi]  -  0 


=  ZZt««xi 


+ 1  = 


The  origin  of  the  coordinate  system  is  selected  at  the  middle 
of  the  beam,  so  that  the  above  three  relations  hold  at  x  equal 
to  plus  and  minus  L/2.  For  a  solution  to  exist  the  resulting 
12  x  12  determinant  of  coefficients  must  be  zero.  All  elements 


of  the  determinant  are  written  in  terms  of  R. ..  The  R. .  terms 
are  trigonometric  or  hyperbolic  functions  of1J  J 


R.*L 

Z  ’ 


x'  1 

R  '2  — 
Kl  2 


S2  «l' 


tj  Rj'8  j  and/or  t^  Rj' 


The  factor  R2L  is  made  a  parameter  F  and  for  a  given  geometric 
cross  section  and  physical  properties  of  the  elastic  and  visco¬ 
elastic  layers  a  variation  of  F  with  hj  can  be  obtained  by  solving 
the  12  x  12  determinant  on  a  digital  computer.  This  will  be  done 
in  the  near  future. 
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Moment 

Real  part  of  M  =  £  £  p^.X^ 

P11  =  t"min0^“  m2no3/^sin  alX  =  chisin  alX 
3/2 

P-J2  =  [m2en0  isin  a]x  =  q-j 2S^ n  alX 

P] 3  =  [ni^a^+tii2no3//2]cos  a-jX  =  q^cos  a-jX 

p14  =  ["m26no3/^cos  alX  =  ^cos  aix 

pi5 =  tmisrm2(si3"3siM2)+ni26(3si2ti‘tl3)]R4i 
+  t'2)]R42 

=  ^15R41+r15R42 

p16  =  [-fn1^+m«{3s  =  2t]-t]3}+m's(s’3-3s|t'2)]R41 

+  [-m1s^'(s]3-3s‘t]2)-m'B(3s’2t^t]3)]R42 

=  q16R41+r16R42 
P17  =  q15R31+r15R32 


P18  =  q16R31+r16R32 


p19  =  [I”1s^-m^(s^3-3s^t^2)+m^(3s^2t^-t^3)]R61 

+  [-m^'+m-Os^t'-t^3)  +  m^e  (s£3-3s't’2)]R62 

=  q19R61  +  r19R62 

pl,10  =  [-n>it^(s^3t'-^3)41n'8(s'3»3s't^2)]R61 

=  ql,10R61+r1t10R62 
Pl,ll  =  q19R51+r19R52 

P1.12  =  ql,10R51+rl,10R52 
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Imag.  part  of  M  =  J  £  a-X.. 

3/? 

ct n  "  [m2Sn  '  ]  sin  a-|X  =  8-| -j sin  a-j X 

“12  =  t"mino"m2no3//2+m2eno3//2^sin  alX  =  e12sin  alX 
a-|3  =  [»m2Bn03X2]cos  a-jX  =  B^cos  a-jX 

=  [m^no^m2(l-B2)no3//2]cos  a-jX  =  b-j^cos  a-jX 

a15  =  tmit]”m2(3sl2tl“tl3)'Hl120^sl3“3sltl2^R41 
+  [m-jS^-m2(s^3-3s|  t^2)-m2B(3sj2t^-tj3)]R^2 

=  615R41+y15R42 
a16  =  q15R41+Y15R42 
“17  =  e15R31+Y15R32 
a18  =  q15R31+Y15R32 

ot-j  g  =  [m1t^-m2(3s^2t^-t^3)-Ki^B(s23-3s^'t22)3R61 

+  [m1s2-m2(s23-3s2t22)-ni2B(3s22t2-t23)]R62 

a19  =  819R61+y19R62 
al,10  =  q19R61+Y19R62 
al,ll  “  619R51+y19R52 

al ,12  =  q19R51+Y19R52 
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I  l  d-.X1  =  Re^l  part  of  y 


Using  a- 


is 

*no  R1 


C  II- 

S  .  0 


*  [— V  +  y  1  sin  a,X  =  t,.  sin  a,X 
nQ  O+e^)  1  11  1 


“eno 

[- - H  Sin  a,X  =  i,«sin  a,X 

(1+B“)  1  12  1 


S  o 

-  o  ]  cos  a  V  =  *  cos  a,X 
nQ  1+^  1  1J  1 


V6 

[ — -]  cos  a  V  =  i  cos  a,X 
1+0^  1  14  ' 


.  r  Ssi  si  et:  st:  t:  es“ 

15 '  + 


=  A15R41+f15R42 


st]  tj  esj  -Ssj  sj 

d16  =  [~T  72  +  r-T  +  -T^4i+[-T  7  +  — 

S1  +t3  1+»  1+g2  4  sj2+tj2  1+0 


=  £16R41+f16R42 


,  ,  Ssi  S1  6ti  ,  stT  t;  0S' 

d17  ”  ^ — ? — 7  ~  — 7  +  "‘ '  — 7t  +  — -  a1  +  — 

s.+ti2  l+a2  i+c^  ^  <•»  2,  ii 2  i .  „2 


Sj  +tj  1+0"  1+06  si  1+6  1+0 


3  *15R31+f15R32 
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> 


d18  =  *16R31+f16R32 


=  *19R61+f19R62 


St  ^ 

dl'10  [s’2+t“2 


3sA  -SsA  sA 

+  Ti7 


'2]R62 


dl,10  =  *1,10R61  +  fl,lOR62 
dl,ll  =  *19R51+fl,SR52 
dl  ,12  =  £l,10R51+fl,10R52 


i 


Imaginary  part  of  y 


Im  y  =  II  9jiXi 


factor  out  —U 

R* 

no^ 

9 1 1  =  [-^y]  sin  a-jX  =  h^sin  a^X 
1+8  , 


g12  =  [; 


Bn, 


1+6 


j]  sin  a-jX  =  sin  a^X 


-ri'a 


9 1 3  =  ( — ^j)  cos  a^X  =  h^3  cos  a^X 
S  no  6 

flu  =  t - V  +  — 9]  cos  a,X  =  h,.  cos  a,X 

14  nQ2  1+82  1  14  1 

Sts  t\  6S '  Ssl  s' 

9 1 5  ~  t  ?  2  ”  7  +  1  ■■■""«"] R«-.+[ . a  a  -  •— -  -a' 

s]2+t]2  1+82  1+82  41  sj^+tp  1+82 


6t' 


1+3 


]R 


42 


=  hl5R4i+ki5R42 

Ss 


3ti 


pi-i 

,J16  =  -  TU  -  7Il'R4l+t 


Sp+tp  1+Bfc  ]n  ’T'  sjfc+t»fc  1+8* 


41 

1  H 


sti  +M 


1+8 


■JR 


42 


=  h16R41+k16k42 


917  "  h15R31+k15R32 
818  =  h16R31+k16R32 


St'  t£  8s; 

9l9  =  +  sf«7 


o  SS« 


c  * 
S2 


34 


1+8  1+8 


OR 


62 
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g19  * 

gl ,10  = 

gi,n  = 

gl ,12  = 


h19R61+  k19R62 


[ 


Ss^ 

Si2+^|2 
s2  z2 


c  1 

s2 


1+8* 


8t*  St' 

- 5'^R6i+f — T — 2 

1+6^  &1  S^+t'2 


4 


1+8  1+8 


2]R62 


hl,10R61+kl,10R62 


h19R51+k19R52 


hl,10R51+kl,10R52 
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Real  part  of  5"  =  J  u^Xj 

uH  =  -nQ,cos  a-|X  =  v-|-|Cos  a-jX 

u12  s  0 

u13  "  -Vin  al^  =  vns1n  ai^ 
u14  '  0 

u15  =  (Si2’r|2^R3r2sltlR32  =  v15R31+w15R32 
u16  =  2sltlR31  +(si2”ti2)R32  =  “  w15R31+v15R32 
u17  =  (sl2”tl2)R41”2sltlR42  n  v15R41+w15R42 
u18  =  2sltlR41+^sl2“ti2^R42  =  ”W15R41+v15R42 
u19  =  {s22“t22)R5r2s2t2R52  =  V19R51+W19R52 
U1 ,10  =  2s^t^R51+(s^2-t^2)R52  =  “wi9R5-|+vi9R52 

u~ll  =  ^s22"t22^R6r2s2t2R62  =  V19R61+W19R62 
U1 ,12  =  2s2t2R61+^s22’t22^R62  =  "W19R61+V19R62 


1 


»  V  .  4 


Imag0  part  of  =  J1..X. 

I J  J 


112  =  -nocos  a-jX 

’is"0 

]14  =  “nosin  alx 
]15  =  w15R31”v15R32 
]16  =  v15R31+w15R32 
^17  w15R41’v15R42 

]18  =  v15R41+w15R42 
]19  =  w19R51"v19R52 
'1,10  =  v19R51+w19R52 
b.ll  =  W19R61"v19R62 
1 1 ,12  =  v19R61+w19R62 


-  29  - 


III.  Determination  of  Natural  Frequencies  -  Approximate  Method 


1.0  In  view  of  the  complexity  in  solving  the  12  by  12  determinants 
for  various  boundary  conditions,  an  approximate  procedure 
was  considered.  The  characteristic  equation  of  the  sixth 
order  homogeneous  differential  equation  developed  in  refer¬ 
ence  5  is  re-examined.  This  is  re-written  as  equation  (III-l). 
When  solved  for  the  natural  frequency,  equation  II 1-4  is  ob¬ 
tained.  An  investigation  of  equation  II 1-4  reveals  that  it 
is  in  the  usual  form  for  homogeneous  beams 


a  4 

2  =  n 


wi  = 


L1* 


tfi-i 


u 

in  which  the  eigenvalue  A0  is  comparable  to  ~y  ,  and  El, 
the  stiffness  in  the  homogeneous  beam,  is  ^  comparable 
to  (B]+B3+K-j62a).  The  values  of  an  are  determined  by  the 
boundary  conditions  and  a  table  of  such  values  is  given 
on  page  35.  The  quantity  (B1+B3+K^62a)  represents  the 
effective  stiffness  on  the  vibrating  beam,  in  which  Bj  and  B3 
are  stiffnesses  of  the  individual  steel  layers  about  their 
own  neutral  axes  and  K-j 62  is  the  portion  of  the  stiffness 
due  to  transferring  the  area  moments  of  inertia  of  the 
elastic  layers  to  the  composite  neutral  axis.  The  factor 
a  is  a  factor  which  indicates  the  shear  carrying  capacity 
of  the  viscoelastic.  When  a  is  equal  to  one-half  then 
there  is  no  shear  strain  in  the  viscoelastic,  whereas  when 
a  is  equal  to  zero,  the  viscoelastic  cannot  transmit  any 
shear  stress  and  each  elastic  layer  bends  independently 
except  for  being  restricted  to  moving  laterally  the  same 
amount.  For  the  case  of  the  simply-supported  beam  having 
no  axial  constraints  at  each  end,  the  value  for  the  eigen¬ 
value  A0  is  exact;  i.e.;  A0  =  (j— )2,  so  that  substituting 
this  value  for  A0  into  the  frequency  equation  yields  the 
exact  natural  frequencies.  The  factor  a  is  exact  for  this 
case  and  truly  represents  the  effect  of  frequency  on  the 
stiffness  of  the  beam.  Guided  by  this  form  of  the  frequency 
equation  and  by  the  exactness  of  using  this  form  for  a  simply- 
supported  beam,  it  is  postulated  that  this  form  of  the  equa¬ 
tion  for  the  natural  frequency  may  be  used  as  an  approxima¬ 
tion  to  the  natural  frequency  of  beams  having  other  end 
conditions.  Thus  it  is  assumed  that  A0  =  bl  where  an  is  the 

usual  factor  determined  by  the  boundary  ^  conditions. 

Thus,  given  the  geometry  of  the  cross-section,  length  of 
the  beam,  the  physical  properties  of  the  materials  and  the 
boundary  conditions,  the  natural  frequency  may  be  approxi- 
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mated  by  equation  I I 1-4  in  which  a  is  given  by  equation 
III-5.  For  ease  of  calculation  the  curves  of  a  versus 
6  are  plotted  for  the  first  five  modes  of  the  cantilever, 
simply-supported,  free-free,  fixed-fixed,  and  fixed-pinned 
beams . 

It  is  recognized  that  this  procedure, suggested  for  finding 
the  natural  frequencies  of  laminated  beams^is  approximate 
and  that  its  accuracy  can  only  be  checked  by  exact  solu¬ 
tions  of  the  kind  performed  in  section  II,  and/or  by  tests 
of  actual  beams.  It  is  felt  that  the  procedure  should 
yield  good  engineering  results. 
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1.1  Approximation  of  Natural  Frequencies 
Equation  (36)  of  reference  3  .is 


fU2d,S+l) 
X4  +  X3  [-L-J - 

0  0  u 


‘1 


]  +  X2  [. 


SR^l+B^O+Sd,) 


h 


1 


1 


d2 

37  “i1 


2SR-,d0u?  d0 

•y-4,— i  -  ^  iW’***)  ■  ° 


III -1 


Solving  for  w^,  we  obtain, 


u2  -  1L  ,2 

W1  c£  Xo 


R-j  i 

>2  +  X^SRj+jj-]  ♦  RfS(l+jS2){fr-+S) 

A2  +  2A  SR,+S2R?(1+B2) 

0  0  I  I 


III -2 


or, 


01 


B,+B,  f  K,«2R1[X  +R-.S  ( 1  +  32)] 
2  =  J_2  A2  1+J  1  0  1 

1  P  o  l 


(B1+B3)[x|  +  2XoSR1+S2R2(1+B2)]  j 


Bl+B3 

where  d,  =  - 

1  k^2 


1 1 1-3 


A  =  _£ 
d2 


K^2 

The  last  equation  may  be  written  as, 


X  B,+B0  K-,62 

a\  =  X2  -=—-=•  +  JL~  a 
1  o  L  p  p 


III -4 


where 


a  - 


Ri[xo+Ris(1+e2)T 


X2  +  2x  SR+S2R?  (1+32) 
o  o  I  i 


1 1 1-5 


If  X  =  n  R, 
o  o  1 


then 


a  = 


n0+S(l+B2) 


n2'+  2nQS+S2{l+32) 


1 1 1-6 


For  a  sandwich  beam  having  a  thin  viscoelastic  layer  S  -  2,  so  that 

-  •  I I 1-7 


.  n0+2(l+62) 

a  =  1 


n0  +  4n0+4(l+32) 
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Reconsidering  equation  HI-5,  for  a,  it  is  seen  that,  for  a  sandwich  beam,  if 
we  let 


then 


or  letting 


anR1l1 2 3+2R12L4 5(l+e2) 

a^+4anR1L2+4R,JL4(1+B2) 


G  =  F2  =  R,L2 

a  G+2G2(l+62) 

d  =  - * - - * - 5J~  . 

an  +4anG+4G  (1+B  ) 


1 1 1-8 


III -9 


III-10 


For  the  special  case  in  which  0  <<1 


1  G 

“  "  T+2 G"  III-ll 

n 


Using  the  relations  III-10  and  III-ll,  plots  of  a  versus  G  are  obtained  and 
are  given  for  the  following  end  conditions 

1.  Cantilever 

2.  Simply-supported 

3.  Free-free 

4.  Fixed-fixed 

5.  Fixed-pi nned 
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For  a  simply  supported  beam  with  no  axial  constraints  on  the  layers  at  each 
support  it  is  found  that 

x0  =  ai =  {r)2  {See  page  35  )  HI-12 


so  that 


(On 


n  4  fG 

-  (i2-)  V 

L  V 


3l+33+^l^a 


III-13 


It  can  be  seen  that  the  term  [B-j+t^+KiS^a]  represents  the  stiffness  of  the 
laminated, beam.  When  is  infinite  (no  shear.strain  in  the  V,E,  layer) 
then  a  -* (or  for  a  sandwich  beam  in  which  S  =  2  then  ct->-  £•)  so  that  the 
stiffness  approaches  a  value  of 


B,  + 


This  is  the  stiffness  of  the  composite  cross-section  in  which  the  shear  is 
carried,  directly  through  the  V.E.  layer  without  shear  deformation.  The 
other  extreme  case  occurs  when  R]  0.  In  this  case  o  becomes  zero,  and 
the  stiffness  is  <JB-|+B3)  and  each  beam  contributes  solely  the  stiffness  about 
its  own  axis;  no  shear  stress  is  'transmitted  through  the  viscoelastic  layer 
although  the  theory  imposed  the  condition  that  both  beams  move  laterally  by 
the  same  amount. 

The  above  description  indicates  that  the  parameter  a  is  an  inverse  measure 
of  the  shear  strain  occurring  in  the  Viscoelastic  layer.  We  note  that  for 
small  values  of  6 


*  1 


so  that  a  maximum  value  of  a  is  ^  for  n0  equal  to  zero;  and  a  approaches  zero 
as  nQ  becomes  large  compared  to  S.  The  quantity  n0  is  always  positive. 

It  is  seen  that  for  the  condition  of  simply-supported,  unrestrained  ends, 
using  \Q  =  (2i.)2,  the  value  obtained  by  the  usual  homogeneous  beam,  an ‘exact 
solution  is  obtained  by  using  equati  on(  1 1 1 -1 3  ).The  factor  a  is  exact  in  this 
case.  This  procedure  suggests  a  means  for  obtaining  the  natural,  frequencies 
laminated  beams  having  other  restraints,  i.e.;  use  the  an  constants  obtained 
for  homogeneous  beams  in  the  relation  X0  =  fg.  ,  A  table  of  values  for  an  is 
given  on  page  35  ,  for  several  boundary  conditions. 


of 
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Beam  Condition 


.  Cantilever 


2.  Simply  - 

supported 


3.  Free  -  Free 


a 

* 

n 

General 

Specific 

Relation 

Relations 

aR  =  (n-^)2Tr2 

a-j  =  3.52 

for  n  >  2 

a2  =  22.0 

a3  =  61.7 

a4  =  121.0 

* 

ag  =  200.0 

i 

j 

an  =  (nn)2 

a-j  =  9.87 

a2  =  39.5 

a3  =  88.9 

a4  =  158. 

a5  =  247. 

aR  *  (n+^)2iT2 

a-j  =  22.0 

a2  =  61.7 

[same  a  for  Free  *• 

a-  =  121.0 

Free  and  Fixed  - 

3 

Fixed] 

a4  =  200.0 

a5  =  298.2 

2  2 

5.  Fixed  -  pinned  an  =  {n+h)  tt 


a-j  =  1504 
a2  =  50,0 
a3  *  104.0 


a4  =  178.0 


a5  =  272. 


Logarithmic,  3  Cycles  X  10  to  the  inch 
Elli  lines  accented. 


G£antmmt,  3  Cyc*c^  X  10  to  the  inch 
5th  lines  accented. 


ogarithmic 


o 


^Hjffwnijrri 


v 

KS  *  ** 


i 


~-*s 
.  i 


urr; 


i 


IV,  Miscellaneous  Considerations 

1,0  The  approximate  expression  for  the  natural  frequency  of  a 
sandwich  beam  further  approximated  for  the  case  in  which  b 
The  expression  for  a  simplifies  to 


2o0 


a  = 


1 


V2 


The  expression  for  cq  is  written  explicitely  in  terms  of  the 
geometric  and  physical  properties  of  the  elastic  and  visco* 
elastic  material* 

Consideration  is  given  to  designing  a  beam  for  optimum  damping 
if  the  physical  properties  of  the  materials  and  the  natural 
frequency  of  the  beam  is  selected*  The  Nation  is  based  on  the 
case  in  which  B  «lo  Considering  <qg  as  the  natural  frequency 
obtained  when  a  =  0,  a  relation,  equation  (V-13},is  obtained 
for  nQO  For  optimum  damping  ref*  5  shows  that  n0  =  5*  Using 
this  tact  a  relation,  equation  (V  =  17),is  obtained  for  the 
necessary  half  thickness,  H-j,  which  would  yield  optimum  damping* 
The  optimum  length  is  obtained  from 


xo  =  noRl  =  jj 


— » 
,  a 


.  i 


yielding 

L  . 

L  ‘•5R-J-*  6 

3,0  Based  on  the  results  of  section  III  ,  for  6  <*1,  expressions 
are  derived  for  G-[  and  3  which  allow  one  to  solve  these  pro¬ 
perties  of  the  viscoelastic  material  if  the  geometric  pro* 
perties,  elastic  properties,  and  the  end-conditions  are  known 
and  if  n  and  <q  are  measured.  These  equations  may  be  used  to 
find  the  shear  properties  of  a  viscoelastic  material  if  b.«K 
This  procedure  may  be  considered  as  a  means  for  testing  material 
to  insure  adherance  to  procurement  specification  of  a  viscoelastic 
material , 

4,0  Using  the  results  of  ref,  5,  an  approximate  expression  for  the 
composite  loss  factor,  equation  (V~23),is  obtained  in  terms  of  h, 
and  a.  It  is  shown  that  this  expression  yields  results  which 
are  within  10%  accuracy. 


f 
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Trefflffa: 


5.0  It  was  shown  in  ref.  3  that  the  natural  frequency  of  a  sand¬ 
wich  beam  increases  as  the  loss  factor,  b,  is.  increased.  It 
is  shown  herein  that  the  same  situation  occurs  using  the 
approximate  expressions  for  the  natural  frequency  developed 
in  section  III,  thus  reinforcing  the  analytical  result  that 
viscoelastic  damping  causes  an  increase  in  natural  frequency 
for  a  laminated  beam. 
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1.1  Detenmi nation  of  the  natural  frequencies  for  small  values  of  material  loss 
7aictor  »  (6  <  1)  . 

Considering  equation  (1 1 1-4 for  small  values  of  0  we  note  in  the  equation 

2  -  x2  rBl+B3  +  V2  ,  TV  , 

“1  "  xo  IV~1 

that  for  a  sandwich  beam  having  a  thin  viscoelastic  layer 

n  +  2(1+b2) 

o  =  -S - k-  ,  IV— 2 

no  +4V4(1+B  > 

and' that  for  0<1 


rT+Z* 

o 


We  note  that  since  n0  can  have  values  from  zero  to  infinity  then  a  has  values 
between  one  half  and  zero.  ' 


Using  the  relation 


an 

_  r)  i  n 
>0  '  noRl  ■  jT 


IV-4 


in  equation  (lV-3) 


then 


.  R^V 


V2RiL 


or  since 


G-j  b  G1b 

R1  =  =  2H2Eb2H1 


IV-6 


then 


4aH9H1E+2G.I 
n  2  1  1 


IV-7 
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-  ’t  . 


or 


a  = 


1 


Now 


- 20 HT  r - 

k,  +  23 


B1  “  TF  )3b  =  |  b  E 


and  for  the  sandwich  beam 


the rf ore 


also 


for 


B1  =  B3 


B1  +  B3  =  |  H^bE. 


K^2  =  Eb2H][2H1]2  =  8EbHB 


H2  «  Hr 


The  value  for  p  is 

Y-j  [  4H-j  b  ]  +  y2[2H2b] 


where 


IV-3 


IV-9 


density  of 


y,  is  the  mass  density  of  the  elastic  material  and  y2  is  the  mass 
the  viscoelastic  material.  The  natural  frequency  may  be  written  as 


a2 

2  _  an 

wi  -  pr 


4HBE 


8eh; 


3t4Y-jH1+^Y2^2^  °  ^yT^T+^y2^2 


IV-10 


or 


2  .  4EH?  4 

Ui  - 


[  l+6ct] 


J1  "  3(4y-|H-|+2y2H2)l4 
in  which  a  is  found  from  equation  { IV~8) . 


IV-11 
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V 

•Sfc  - 


Thus,  in  order  to  calculate  the  natural  frequency  of  a  sandwich  beam,  one 
would  need  to  know 

1.  a  -  This  is  determined  from  the  end  conditions  and  the  mode  of 

n  vibration.  A  list  of  these  values  is  given  on  page 35. 

2.  (2H-|)  -  This  is  the  thickness  of  the  elastic  layers. 

3.  (2H2)  -  This  is  the  thickness  of  the  viscoelastic  layer. 

4.  L  -  This  is  the  length  of  the  beam. 

5.  E  -  Young's  Modulus  of  the  elastic  material. 

6.  G,  -  This  is  the  real  part  of  the  shear  modulus  (storage  modulus) 

of  the  viscoelastic  material. 

7.  y-|  -  This  is  the  mass  density  of  the  elastic  material. 

8.  y?  -  This  is  the  mass  density  of  the  viscoelastic  material 

L  (Note  that  the  mass  per  unit  length  and  per  unit  width  is  the 
quantity  4y1H1  +  2y2H2). 
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2 j  Design* ~of -a  Sandwich  Beam  for  Optimum  Composite  -  Loss  Factor 

Using  the  suggested  approximate  procedure  for  finding  .the- natural  fre¬ 
quency  of  a  sandwich  beam  having  a  thin  viscoelastic  layer,  we  have 

•  an 

xo  =  noRl  *  • 

But  using  eq.  (  ). 


?  a  2  4E  H,3 
2  _  n  r«  I 

uio "  r?  [- 


*3 


31*  4y  -j  H-j  +  2y2H2 

we  can  obtain  from  ( IV- 12) 

,  „  *2  ’  3(4ylHl  +  2y2H2^  “10 

(n  r  )  - - - - 

01  4  EH? 


IV- 1 2 


IV-13 


For  optimum  damping, 

h-j  =  50  (see  pg.  11-33  ref.  5  ). 

Also  for  h1  =  50  we  find  from  the  nQ  vs  h]  plot  that  nQ  =  5. 


for  nQ  =  5,  we  find  for  3«1  that 
2  2  r13n 

W1  *  w10  [T] 


or 


2  7  2 

wio '  TJ  “r 


IV-14 


IV-15 


Since  Ri = 

then  9 

G1 

25[  j — j]  = 

16E^ 


21  (4y1H1  +  2y2H2)^ 


52  EH 


IV-16 
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Solving  for  H-j  we  find 


H 


^2^2^wl 

2  2  2  ' 
3,86  -  4y  1 


IV-17 


Thus  for  given  elastic  material  (E  and  y, )  and  a  given  viscoelastic 
material  (G,  and  y2)  one  may  select  an  w,  witn  its  associated  G-j  and  a 
viscoelastic  thickness  layer  (2H2),  then  solve  for  the  half  thickness  of 
the  elastic  layer  (H-j)  for  optimum  damping,,  The  length  of  the  beam  may  be 
found  from  the  expression. 


%Ri  * ;? 


so  that  for  nQ  =  5 


n  -h 


L  = 


4a  EH,H,  . 

r  n  1  2n% 

[ — 5^j  3 


rv-18 


in  which  an  is  found  on  page  35. 

Equations  (IV-17)and (lV-18)yield  the  half  thickness  of  the  elastic  layer 
and  the  associated  length  of  trie  beam  to  achieve  optimum  damping. 


I 
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3.1  A  Suggested  Method  for  Determining  the  Shear  Moduli  and  Loss-Factor  of  a 
r  Vi  scoelas ti  c-  MaterfaT  -  B«  C 


A  method  for  determining  the  shear  moduli  and  loss-factor  of  a  visco¬ 
elastic  material  is  suggested,  based  on  the  analysis  on  pages 44  and  45 
If  in  some  manner,  as  for  example  impedance  measurements  of  a  free-free  beam, 
the  composite  loss  factor  n  and  the  natural  frequency  <01  and  its  associated 
mode  number  are  measured  for  a  sandwich  beam  having  a  thin  viscoelastic  layer, 
then  the  shear  modulus  and  material  loss  factor  8  may  be  calculated,, 

For  the  case  in  which  8  is  small  G-j  and  n  can  be  obtained  in  the  follow¬ 
ing  manner.  Knowing  the  mode  number  and  the  end-conditions,  the  value  of  an 
may  be  obtained  from  table  1„  This  may  be  used  in  equation  (lV-110-to  solve 
for  a,  i,  e,, 


a 


3(4y]H^2y2H2)L4  2  1 
24E  a^  ^  6 


IV-19 


Using  equation  ( IV-7)and  solving  for  G-j ,  one  obtains 


G 


1  \J~[  1  —  2a  ] 


t 


IV— 20 


in  which  a  is  obtained  by  solving  eq.  ( IV-.l 9) .The  real  part  of  the  shear 
modulus  is  thus  calculated. 


The  relation  for  the  composite  loss  factor  n  for  a  sandwich  beam  is  shown 
in  ref.  5  ,  to  be  eq.  11-15. 


¥1*0 

((R,s+x0)2+(r1sb)2] 


IV-21 


which  for  the  assumptions  being  considered  i.e.  (S<1),  becomes, 
(4y1H1+2y2H2)  »i[G1+2X0  H2H,E]2 

3  =  ri  - - - - - 


IV-22 


This  last  equation  allows  one  to  calculate  the  material  loss  factor  8. 
Having  found  G-j  and  8,  the  loss  modulus  G2  may  be  found  using 


G 


2 


8G 


r 
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v. 


V  4,1  An  Approximate  Equation  for -the  -Composite  Loss-Factor 

It  was  found  in  ref.  5  ,  that  for  a  sandwich- beam  having  a  thin 
viscoelastic  layer 


6  "o3* 

h,[(2+n0)2+(2s);!] 


An  inspection  of  the  n  vs  lr^  plot  for  values  of  6  between  .1  and  1 
(pg '11-34  ref.  5  )  indicates  that  the  log  n  vs  log  h-j ,  curve  is  approximately 
a  straight  line  and  rather  independent  of  6  for  values  of  $  between  .1  and  1. 
The  relation  can  be  assumed  to  be 


n 


0 


=  c 


IV-23 


in  which  c  and  k  must  be  evaluated  from  the  curve*  We  see  that  two  corres¬ 
ponding  points  on  the  nQ-h-|  plot  are 

n  r  20  and  h-,  =  500 
0  l 

n  =  2  and  h,  =  11 
0  1 

Solving  for  c  .and  k  by  using  the  above  . values  we  find 

k  =  06 

and  c  =  .475 


This  then  yields  the  relation 


n  =  .668 


[(2+.475h1c6)2+(2B)2] 


IV-24 


A  check  of  three  values  of  h]  indicates  a  good  correlation  with  the  exact  re¬ 
sults.  These  are  shown  below. 


n  __  I 

hl 

Approx. 

Exact 

%  error 

8 

4 

.022 

.0238 

-7.5 

.1  • 

4 

.153 

.150 

+1.95 

1  ' 

40 

.031 

.033 

-6.05 

.1 

40 

.282 

.280 

+  .7 

1 

600 

.0195 

.0185 

+5.4 

.1 

600 

.195 

.185 

+5.4 

i 

1 

) 
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The  approx,  value  of  n  uses  equation  ( IV-»24)  whereas  the  exact  value  of 
n  is  taken  fron?  the  n  vs  h-j  plot  of  ref.  5.. 


As  was  shown  previously  in  ref.  5  ,  the  approximate  equations  for  small 
values  of  h-j  i.e.  hi  less  than  1,  and  the  approximate  equation  for  larger  values 
of  h^t  i.e.,  h-j  >  2000  are  given  as 


.  3h^ 

n  =  ~nr 


for  h-j  <  1 


n 


for  h-|  »  I. 
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5j  Comments  on -the  apparent  increase  in  w-j  when  damping  is  present 


In  refo  2  ,  it  was  shown  that 


Ui  -  0), 


[1  +  ¥•] 


where 


w-j  is  the  natural  frequency  including  damping  effects 
2  EI^ 

uQ  =  p  [ — ] 2,  the  undamped  natural  frequency 
n  is  the  composite  loss-factor 
p  is  associated  with  the  mode  number 


El  -  stiffness 
p  -  mass/unit  length. 

In  the  above  El  is  the  stiffness  and  should  take  into  account  the  decrease 
in  stiffness  due  to  the  shearing  effect  occuring  in  the  viscoelastic  layer., 

)  Thus  if  used  properly  I  is  the  total  I  of. the  cross-section  about  the  composite 
neutral  axis,  at-the  low  frequencies,  and  decreases  to  the  sum  of  the  individual 
I's  about  their  own -neutral  axis  as -the -vibration  frequency  increases.  But  for 
a  given  El,  indeed  the  relation  above  indicates  the  natural  frequency  would 
increase  with  the  addition  of  damping.  As  a  further  substantiation  of  this 
increase  we  look. at  the  expression,  for  the  sandwich  beam  having  a  thin  visco¬ 
elastic  layer, 

U1  =  “lo  ^  +  6<x^ 

The  factor  a  accounts  for  the  shear  effect  of  the  viscoelastic  layer  as 
was- pointed  out  previously  since  ,  for  0-0,  it  varies  from  0  (no  shear 

carrying. capacity  by  V«E,)  to  h  (all  shear  carried).  Now  a  also  contains  3, 
the  material  loss  factor.  In  particular  we  have 

nQ  +  2(1  +  82) 

y  -  q  p 

n  +  4nQ  +  4(1  +  Bc) 


If  we  compare  a  containing  0  and  that  for  which  0  =  0  i„e„  cx  ,  ana  it  cms 
ratio  -§•  is  greater  than  one,  then  the  natural  frequency  of  a  beam  would  tend 
to  be  higher  with  damping  than  it  would  be  without  damping.  Looking  at  this 
ratio  we  see,  n  +  2  (i 


and  if  this 


n. 


+  e2) 


•£L 

ct. 


no2  +  4nQ  +  4TV'+  02)  (no+2)[no+2(l+02)] 


n  +T~ 


n^+4  n0+4  (1+02) 
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2  n  B 


n0  +  4  no+40+fl4‘) 


For  e  small  compared  to  nQ,  we  have 


2-  =  1  + 


O  O 

in  +  2)2  B 


This. shows  2->  therefore- this  indicates  that  the  natural  frequency  will 
ao 

increase' with- damping,  from  the  value  it  would  have  with  no  damping. 
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V.'  Errata 

The  following  two  equations  are  the  errata  for  equations  H-1 
(Moment  and  shear)  of  reference  5. 

Let 

jm-j3Ki(S2+(Bi+B3)S 

II'9  B,+B, 
jm9* - 7- 

V  1  ^(i+r) 

Then, 

3 

(11-10)  *  E^[-m^a^sin  -m2ai  sin  a^x] 

+  Ei2^m26a13s *n  alx^“E22tm2pa13cOS  a1X^ 

+  E21[m-ja.jCOS  a^x  +  m2a^  cos  a-jX] 

+  E31tmlS31“m2U3l'Hn2BUi21+E32["‘nlS3<:+m2U32+m2eU3T;l 
+  E41[m1S41-m2U41-Hn2BU42)+E42[-m1S42+in2U42+<n2BU41] 

+  E51[m1S5rm2U51-*m2BU52]tE52[-m1S52-hn2U52+m2BU51] 

+  E61(«J1S51-m2U614«>289621+E62["nlS624ffl2U62',,n26U611 

+i^1-1[m2eai3sir!  a1x)-E21[m2ea13cos  a^] 

+  E12;-m1a1  sin  a^x-fn^sin  a^+m^a^sin  a^] 

3 

*  E22tm1a1cos  a^-hrija^cos  a,x-m2sa,  cos  a,xj 
+  E31  (»1  S32-"'2932'm2BU31 ] +E32tml  S3l  -m2U31+m28U321 

+  E41[m1S42-m2U42+ni2BU41)+E42(in1S41-m2U41«2BU42] 

+  E51(m1S52-m2U52+m2SU51]+E52(m1S51-ni2U51+m2sU52) 

+  E61["11S62'™2U62+n'2BU61)«62imlS6rm2U6l'iTO26U625] 


and  11-11 


11-15 
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Evaluation  of  Shear  -  V 


r 


2  4 

(11-11)  -  $•*  cosa^x-m2ai  cosa^x] 

4  4 

+  E12[m2ea.j  ,cps  aixl+E22^m2eal  sin  a-jX] 

2  4 

+  Eg^t-m^a^  sin  a^x-m2a-|  sin  a-jx] 

+  ^3^  tmiT3i“m2V31+m2BV32^+E32^"mTT32‘4'm2V32'fni2BV31 

+  E4*|  tm*|'r4-j~m2V41+m20V423+E42E-m^T42-Hn2V42-Hn23V41 3 

+  E51^mlT5rm2V51+m26V52^+E52^mlT52+m2V52+m2eV51^ 

+  E61^mlT6rrn2V614in2^V62^+E62^“mlT624tI,2V62'Hn23V61^ 

+^E11[m28a<14cos  a-jX]+E2^[m2$a-j4sin  a^x] 

+  E^t-m^a^cos  a-jX  -  m2a-j4cos  a^x  -Hn^a^cos  a-jX] 
+  E22[-m^a-j2sin  a-jX  -  m2a^4sin  a^x  +rn26a4sin  a-jx] 

+  E3^C>niT32-m2V32+m2BV31l+E32[m1S3^-m2V314tn2BV32] 

+  E41  tmiT42*'m2V42'HT12BV41  ^+E42^ml  S41  ”*m2V41+m2BV42-* 

+  Eg-jtm^Tgg-mgVgg+ingBVg^+Eggtm^g^-mgVg-j^BVgg] 

+  E6ltmi’r62"ni2V62+m28V61^+E62^mlS6rnl2V61+ni2eV62'J?' 


11-16 

t 
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